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ABSTRACT 

The spectral and temporal behavior of exoplanet host stars is a critical input to models of the 
chemistry and evolution of planetary atmospheres. Ultraviolet photons influence the atmospheric 
temperature profiles and production of potential biomarkers on Earth-like planets around these stars. 
At present, little observational or theoretical basis exists for understanding the ultraviolet spectra of 
M dwarfs, despite their critical importance to predicting and interpreting the spectra of potentially 
habitable planets as they are obtained in the coming decades. Using observations from the Hubble 
Space Telescope^ we present a study of the UV radiation fields around nearby M dwarf planet hosts that 
covers both FUV and NUV wavelengths. The combined FUV-I-NUV spectra are publically available 
in machine-readable format. We find that all six exoplanet host stars in our sample (GJ 581, GJ 
876, GJ 436, GJ 832, GJ 667C, and GJ 1214) exhibit some level of chromospheric and transition 
region UV emission. No "UV quiet" M dwarfs are observed. The bright stellar Lya emission lines are 
reconstructed, and we find that the Lya line fluxes comprise ~ 37 - 75% of the total 1150 - 3100 A flux 
from most M dwarfs; > lO'^ times the solar value. We develop an empirical scaling relation between 
Lya and Mg II emission, to be used when interstellar H I attenuation precludes the direct observation 
of Lya. The intrinsic unreddened flux ratio is F(Lya)/F(Mg II) = 10 ± 3. The F(FUV)/i^(NUV) 
flux ratio, a driver for abiotic production of the suggested biomarkers O2 and O3, is shown to be 
~ 0.5 - 3 for all M dwarfs in our sample, > 10'^ times the solar ratio. For the four stars with moderate 
signal-to-noise COS time-resolved spectra, we find UV emission line variability with amplitudes of 
50 - 500% on 10^ - 10^ s timescales. This effect should be taken into account in future UV transiting 
planet studies, including searches for O3 on Earth-like planets. Finally, we observe relatively bright 
H2 fiuorescent emission from four of the M dwarf exoplanetary systems (GJ 581, GJ 876, GJ 436, 
and GJ 832). Additional modeling work is needed to differentiate between a stellar photospheric or 
possible exoplanetary origin for the hot (T(H2) ~ 2000 - 4000 K) molecular gas observed in these 
objects. 

Subject headings: planetary systems — stars: individual (GJ 581, GJ 876, GJ 436, GJ 832, GJ 667C, 
GJ 1214) — ultraviolet: stars — stars: activity — stars: low-mass 



1. INTRODUCTION 

Stellar ultraviolet (UV) photons play a crucial role in 
the dynamics and chemistry of all types of planetary 
atmospheres. Extreme-UV (EUV; 200 < A < 911 A) 
photons capable of ionizing hydrogen heat the atmo- 
sphere and are the primary drivers of atmospheric es- 
cape in both the Solar System and exoplanetary sys- 
tems. For short-period gas giant planets, EUV irra- 
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diati on can heat the atmosp here to lO'' K (lYelld 
l2004t iMurrav-Clav et al.ll2009D resul ting in atmospheric 
mass- loss ( Sanz-Forcada et al.l 12011) and the inflated 
exospheres observed in H I Lya (iV i dal-Madiar et al.l 
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120031 : iLecavelier des Etanes et aP I2012D a r id metal ab- 
sorption line fVidal-M adiar et al.l 1200^ iLinskv et al.l 
I2OIO; Fossati ct al. 2010) transit spectra. For terres- 
trial atmospheres, an increase in the EUV flux can el- 
evate th e temperature o f the thermosphere by a factor 
of > 10 (jTian et al.ll200l . potentially causing significant 
and rapid atmospheric mass-loss. 

Far-UV (FUV; 912 < A < 1700 A) and near-UV (NUV; 
1700 < A < 4000 A) photons from the host star are 
the drivers of atmospheric chemistry through their ef- 
fects on the photoexcitation and photodissociation rates 
of many abundant molecular species. Solar chromo- 
spheric emission (specifically H I Ly/3) excites H2 fluo- 
rescence in gas giant planets, producing rovibrationally 
excited H2 that catalyzes chemistry in t he Jovian atmo- 
sphere (|Cravenslll987t IKim fc Foxiil994 ). In exo- Jovian 
atmospheres heated to T > lOOOK, Lya fluorescence of 
both H2 an d CO may become an important excitation 
mechanism (jWolven et al.|[l997l ): this process may be an 
observable diagnostic f or ga s giant atmospheres outside 
the solar system (jYelld 120041: [France et al.ll2010b( ). H2O, 
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CH4, and CO2 are sensitive to FUV radiation, in par- 
ticular the bright H I Lya hne which has spectral co- 
incidences with these species. UV radiation can be im- 
port ant for the long-terrn habitability of Earth-like plan- 
ets ("Bucci no et al.ll2006[) , and the combination of FUV 
and NUV photons can influence the molecular oxygen 
chemistry. For example, the O2 abundance can be signif- 
icantly enhanced through the photodissociation of CO2 
and H2O. The subsequent formation of O3 depends sensi- 
tively on the spectral and temporal behavior of th e FUV 
and NUV radiation fields of the host star (Tia n et"al] 
I2OI2I: iDomagal-Goldm an et al.ll2012l) . 

The UV (FUV -f NUV) spectra of solar-type stars have 
been studied extensively, both observationally and theo- 
retically (e.g.. Woods et al. 2009; Fontenla et al. 2011; 
Linsky et al. 2012c). By studying solar analogs of var- 
ious ages, we also have an und erstandi ng of the evolu- 
tion of G star spec tra (|Avresl [1997; .Riljas et all 120051: 
iLinskv et al.ll2012Hl . By contrast, the observational and 
theoretical literature on the UV spectra of M dwarfs is 
sparse. The photospheric UV continuum of M dwarfs is 
very low relative to solar-type stars due to their lower 
effective temperature. Bright chromospheric and transi- 
tion region emission lines dominate the UV spectrum of 
M dwarfs, comprising a much larger fraction of the stel- 
lar luminosity than for solar-type stars. Outside of a few 
well-studied flare stars (e.g., AU Mic, AD Leo, EV Lac, 
Proxima Cen), M dwarfs have largely been ignored by 
UV observers because most investigations were aimed at 
understanding the origin and nature of energetic events 
on low- mass stars. There are very few UV observations 
of the older, quiescent M dwarfs that are most likely to 
harbor potentially habitable planets, and at present no 
theoretical atmosphere models exist that self-consistently 
predict the UV emission from an M dwarf. 

UV variability of M dwarf exoplanet host stars 
is also essentially unconstrained by observation at 
present. M dwarfs display variability on many 
timescales and flare activity occurs with greater fre- 
quency and a mplitude than on main-sequen ce G 

and K dwarfs (Hawle v et all [19961 IWest et al.l [20P: 
IWelsh et al...2007 : Walkowicz et al.ll201lD . Because most 
flare activity is thought to be related to magnetic en- 
ergy depositio n in the corona, transition region, and 
chromosphere (jHaisch et al.lll99l[) , UV flare activity has 
an amplitude comparable to or greater than that ob- 
served in optical light curves. UV flare behavior on 
active M dwarfs is characterized by strong blue/ NUV 
continuum emission and the enhancement of chromo- 
spheric and transition region emission lines (e.g.. Haw- 
ley & Pettersen 1991; Hawley et al. 2003). FUV flares 
may also be signpost for enhanced soft X-ray, EUV, 
and energetic particle arrival rates on planets orbiting 
these stars, all of which can play an important role 
for the chemistry and evo lution of exoplanetary atmo- 
spheres (jSegura et a l."2010') and the long -term habitabil- 
ity of these worlds (Buccino et al.|[2007lf . 

Given the paucity of existing observational data and 
lack of theoretical models of the UV spectra of M dwarfs, 
the most reliable incident radiation field for photochemi- 
cal models of extrasolar planetary atmospheres is one cre- 
ated from direct observations of the host star. Towards 
the goal of a comprehensive library of UV radiation fields 
for use as inputs in models of exoplanet atmospheres, we 
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Fig. 1. — The observed 1150 - 3140 A fluxes from the M dwarf 
exoplanet host stars in the MUSCLES sample with complete (Lya 
+ FUV + NUV) datasets. The 1760 - 2100 A emission has been 
omitted because the flux levels are below the STIS G230L detection 
level. For display purposes, the data have been convolved with a 
3 A FWHM Gaussian kernel and offset as follows: GJ 581, no 
offset; GJ 876, logioF^ -I- 1.2; GJ 436, logioF^ -|- 2.4; GJ 832, 
logioF^ + 3.6. 

have carried out a pilot program to observe the spectral 
and temporal behavior of M dwarf exoplanet host stars. 
Measurements of the Ultraviolet Spectral Characteristics 
of Low-mass Exoplanet host Stars (MUSCLES). 

In this paper, we present observational results from the 
MUSCLES pilot program, spectrally and temporally re- 
solved UV data for six M dwarf exoplanet host stars (GJ 
581, GJ 876, GJ 436, GJ 832, GJ 667C, and GJ 1214) ob- 
served with the Hubble Space Telescope-Cosmic Origins 
Spectrograph (COS) and Space Telescope Imaging Spec- 
trograph (STIS). We describe the targets and their plan- 
etary systems in Section 2. Sections 3 and 4 describe the 
observations, the spectral reconstruction of the impor- 
tant stellar Lya emission line, and the creation of light 
curves for chromospheric and transition region emission 
lines. In Sections 5 and 6, we describe the spectral and 
temporal characteristics, respectively, of our target sam- 
ple. Section 7 places these results in context of models of 
terrestrial planetary atmospheres and discusses implica- 
tions for UV transit studies. We summarize our results 
in Section 8. 

2. TARGETS 

In this section, we briefly summarize the stellar and 
planetary system properties of our M dwarf targets. Op- 
tical spectra of most of our targets have been published 
in the literature, and while we do not make a detailed 
comparison of the UV and visible emission lines here, it 
is valuable to place the MUSCLES targets in the context 
of traditional optical activity indicators. The stars in our 
sample would traditionally be considered "opticall y inac- 
tive" , based on their Ha absorption spectra (Gizis et al.l 
l2002f ). However, all of our stars with measured Ca II H 
and K profiles show weak but detectable emission (equiv- 
alent widths, EW(Ca II) > 0), indicating that at least a 
low level of chromospheric activity is present in these 
stars ((Rauscher fc Marcv. ,2006; Walkowicz fc HawlevI 
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Fig. 2. — ('e/t) H ST-STIS E140M and G140M Lyo profiles {black histograms) . The two-component Gaussian emission line + interstellar 
absorption model is overplotted in pink. The neutral hydrogen and deuterium absorption components are marked with blue arrows. The 
spectra are scaled as marked and offset in +10 X lO"'^* erg cm~^ s~^ steps, photometric uncertainties are displayed as green error 

bars, (right) Intrinsic two-component Gaussian model to the stellar Lya emission line, scaled to the flux at 0.16 AU, in the habitable zone 
of these M dwarf exoplanet host stars. The spectra are scaled as marked and offset in +20 erg cm~^ steps. 



|2009() . Our stars show Ha in absorption with equiv- 
alent widths in the range -0.4 < EW(Ha) < -0.2 A 
and Ca II emission with equivalent widths in the range 
0.2 < EW(Ca II) < 0.5 A. Adopting the M dwarf clas- 
sification from IWalkowicz fc HawlevT (|2OO90 . these stars 
would be referred to as possessing intermediate chromo- 
spheres, or weakly active M dwarfs. 

GJ 581 ~ GJ 581 is an M2.5 dwarf at a distance 
of 6. 3 pc. It is estim ated to have an age of 8 ± 1 
Gyr (jSelsis et al.ll2007[) and a so mewhat subsolar metal- 
licity, [Fe/H] = -O-IO - -0.02 ([Johnson fc AppsI [200l 
iRoias-Avala et all2010[) . GJ 581 is not detected in X-ray 
surveys (logio^x < 26.89 erg s~^; Poppenhaeger et al. 
2010) and its optical spectrum displays Ha in absorption, 
therefore chromospheric and coronal activity are thought 
to be low for this target. GJ 581 has one of the richest 
known planetary systems, with possibly up to six planets 
(four confirmed) includin g several of E arth/super-Earth 
mass (IM avor et ail [20091 iTuomil [20TlD . GJ 581d is a 
super-Earth (Mp « 6 Mq) that resides on the outer edge 
of the habitable zone (HZ; ap = 0.22 AU; Wordsworth 
et al. 2011; von Braun et al. 2011). 



GJ 876 ~ GJ 876 is an M4 dwarf at 4.7 pc, and is the 
only planet-hosting M dwarf with a well-characterized 
UV spectrum prior to the present work ([Walkowicz et al.l 
[20081: iFrancc et al.ll2012aD . GJ 876 has super-solar metal- 
licity ([Fc/H] = 0.37 - 0.43; Johnson & Apps 2009; 
Rojas-Ayala et al. 2010), and differing estimates on 
the stellar rotation period (40 < < 97 days) result 
in large uncertai nties in the age estimate for this sys- 
tem, 0.1 - 5 Gyr ([Rivera et al.|[2005l [2010l: IGorreia et all 
I2010D . While the star would be characterized as weakly 
active based on its Ha absorption spectrum, UV and 
X-ray observations have shown the presence of an ac- 
tive upper atmosphere (Walkowicz et al. 2008; France 
et al. 2012a; logio^x = 26.48 erg s~^; Poppenhaeger 
et al. 2010). GJ 876 has a rich planetary system, 
with four planets ranging from a super-Earth (GJ 876d, 
Mp w 6.6 M^) in a short-period orbit (ap = 0.02 AU; 
Rivera et al. 2010) to two Jovian-mass planets in the HZ 
(GJ 876b, Mp « 2.27 Mj„p, ap = 0.21 AU; GJ 876c, 
Mp w 0.72 Mj„p, ap = 0.13 AU; Rivera et al. 2010). 

GJ 436 - GJ 436 is an M3 dwarf star located at a 
distance of 10.3 pc. It has a 45 day rotation period, a rel- 
atively old age (~ 6^5 Gyr; Torres 2007), and may have 
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Fig. 3.— C IV AA1548,1550 A resonance doublets for the four 
stars with Lya detections and FUV spectra, plotted in the helio- 
centric velocity frame. The emission lines have been scaled to the 
flux at 0.16 AU, in the habitable zone of these M dwarf exoplanet 
host stars, and offset in 0.3 erg cm~'^ steps. The dashed 

lines are the rest velocities of the 1548 A and 1550 A lines, respec- 
tively. 



TABLE 1 
MUSCLES HST Observations. 



Target 


Date 


Instrument, Mode 


PID 


Texp (s) 


GJ 581 
GJ 581 
GJ 581 
GJ 581 


2012 Apr 18 
2012 Apr 18 
2011 Jul 20 
2011 Jul 20 


STIS, G140M 
STIS, G230L 
COS, G130M 
COS, G160M 


12034 
12034 
12034 
12034 


1107 
265 
1478 
900 


GJ 876 
GJ 876 
GJ 876 
GJ 876 


2011 Nov 12 

2011 Nov 12 

2012 Jan 05 
2012 Jan 05 


STIS, G140M 
STIS, G230L 
COS, G130M 
COS, G160M 


12464 
12464 
12464 
12464 


1138 
267 
2017 
2779 


GJ 436 
GJ 436 
GJ 436 
GJ 436 


2010 Jan 05 
2012 May 10 
2012 Jun 23 
2012 Jun 23 


STIS, G140M 
STIS, G230L 
COS, G130M 
COS, G160M 


11817 
12464 
12464 
12464 


1762 
1665 
3372 
4413 


GJ 832 
GJ 832 
GJ 832 
GJ 832 
GJ 832 


2011 Jun 09 

2012 Jun 10 
2012 Apr 10 
2012 Jul 28 
2012 Jul 28 


STIS, E140M 
STIS, E230H 
STIS, G230L 
COS, G130M 
COS, G160M 


12035 
12035 
12464 
12464 
12464 


2572 
3135 
917 
2163 
2925 


GJ 667C 
GJ 667C 


2011 Sept 04 
2011 Sept 04 


STIS, E140M 
STIS, E230H 


12035 
12035 


2396 
3023 


GJ 1214 
GJ 1214 
GJ 1214 
GJ 1214 


2011 Apr 27 

2012 Aug 12 
2012 Aug 04 
2012 Aug 04 


STIS, G140M 
STIS, G230L 
COS, G130M 
COS, G160M 


12165 
12464 
12464 
12464 


7282 
1620 
3289 
4368 



a super-solar metallicity ([Fe/H] ~ 0.00 - 0.25; Johnson 
& Apps 2009; Rojas-Ayala et al. 2010). GJ 436 does 
show signs of an active corona (logio^x = 27.16 erg s~^; 
Poppenhaeger et al. 2010), an d its chromosphe r ic Lyg 
emission has been observed by lEhrenreich et al.l ()20II|) . 
GJ 436 is no table for its well-studied transiting N eptune 
mass planet (|Butler et al.ll2004l : iPont et al.ll2009f ). orbit- 
ing at a semi-major axis of w 0.03 AU, interior to its 
HZ (0.16 - 0.31 AU; von Braun et al. 2012). Addi- 
tiona l low-mass planets ma y also be present in this sys- 
tem (jStevenson et al.ll20I2D . 

GJ 832 - GJ 832 is an Ml dwarf at d = 4.9 pc. GJ 
832 is not as well characterized as other targets in our 
sample; an age determination for this star is not avail- 
able. Coronal X-rays have been detected from GJ 832, 
logioijf = 26.77 erg s~^ (Poppenhaeger et al. 2010). 
This subsolar metallicity star ([Fe/H] — —0.12; Johnson 
& Apps 2009) hosts a 0.64 Mj^p mass planet in a 9.4 
year orbit (ap = 3.4 AU; Bailey et al. 2009). 

GJ 667C ~ GJ 667C (M1.5V) is a member of a 
triple star system (GJ 667AB is a K3V + K5V bi- 
nary) at a distance of 6.9 pc. Th is 2-10 Gyr 
M dwarf ()Anglada-Escude et al.l I20I2D is metal-poor 
([Fc/H] = -0.59 ± 0.10, based on an analysis of GJ 
667AB; Perrin et al. 1988) and may host as many as 
three planets, including a super-Earth mass planet (GJ 
667Cc, Mp w 4.5 M®, ap = 0.12 AU) orbiting in the HZ 
(0.11 - 0.23 AU; Anglada-Escude et al. 2012). 

GJ 1214 - GJ 1214 is a late M dwarf (M6V) at 13 
pc, making it the coolest and most distant of the targets 
in th e MUSCLES pilot study. It has an age of 6 ± 3 
Gyr ([Charbonneau et all I2009() . a super-solar metallic- 
ity ([Fe/H] = +0.39 ± 0.15; Berta et al. 2011; Rojas- 
Ayala et al. 2010), and show s signs of optical flare ac- 
tivity (jKundurthv et al.|[20T]l ). GJ 1214b is a transiting 
super-Earth (Mp « 6.5 M®, ap = 0.014 AU; Charbon- 
neau et al. 2 009), possibly harboring a dense, wa ter-rich 
atmosphere (jBean et al.l 120 lot IDesert et al.ll201lD . 

3. OBSERVATIONS 

The MUSCLES observing plan uses the two pri- 
mary ultraviolet spectrographs on HST to create quasi- 
continuous M dwarf exoplanet host star data from 1150 - 
3140 A. Data for the MUSCLES program were obtained 
as part of HST GTO and GO programs 12034, 12035, 
and 12464, acquired between 2011 June and 2012 Au- 
gust. COS has a factor of > 10 x more effective area than 
the medium resolution modes of STIS in the FUV, and 
its low detector background and grating scatter make it 
factors of ~ 50 x more efficient for the study of faint far- 
UV chromospheric and transition region emission from M 
dwarfs (Green ct al. 2012). We used COS, in its time-tag 
(TTAG) acquisition mode, to observe the FUV spectra of 
our target stars and study time variability on scales from 
minutes to hours (§4.2 and §6). Multiple central wave- 
lengths and focal-plane offset (FP-POSs) settings with 
the COS G130M and G160M modes were used to create 
a continuous FUV spectrum from 1145 - 1795 A. These 
modes provide a point-source resolution of Av » 17 km 
s~^ w ith 7 pixels per resolution element (jOsterman et al.l 

[gonl). 

Because COS is a slitless spectrograph, observations 
at H I Lya are heavily contaminated by geocoronal 
emission. We therefore observed the stellar Lya pro- 
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TABLE 2 

Stellar Lyo Emission Line Widths and Interstellar, Atomic Hydrogen^. 



Target 


Fjiarrow 1 ^broad 


FWHMnarroTjj 


FWUM^road 


logio N{m) 


bni 






(km s-1) 


(km s-i) 


(cm^^) 


(km s^-*-) 


GJ 581'' 


8.4 


110 


228 


18.35 ± 0.06 


10.1 


GJ 876*^ 


5.0 


132 


303 


18.06 ± 0.03 


8.0 


GJ 436*^ 


11.8 


109 


294 


18.19 ± 0.04 


7.4 


GJ 832<= 


7.6 


96 


163 


18.47 ± 0.02 


9.7 


GJ 667C= 


8.9 


85 


233 


18.07 ± 0.03 


12.7 


AY) Leo*^ 


2.9 


166 


412 


18.47 ± 0.01 


9.0 


HD 189733''''' 


3.5 


170 


480 


18.45 


10 



^ - Lya profiles reconstructed using the iterative least-squares technique described in §4.1. 
b STIS G140M observations. 
STIS E140M observations. 

HD 189733 reconstruction failed to converge, the presented profile parameters provide a fit to the observed line 
profile that is qualitatively similar to the formal fits for the M dwarf spectra. 



file with STIS, using either the G140M/cenwave 1222 
mode {Av w 30 km s^^) through the 52" x 0.1" slit, 
or the E140M mode (Ai; « 7.5 km s"^) through the 
0.2" X 0.2" slit for the brightest targets. In order to in- 
clude the important Lya emission line for GJ 436 and 
GJ 1214, we downloaded STIS G140M observations of 
these stars from the MAST archive (program IDs 11817 
and 12165, respectivel y). The Lya p rofile of GJ 436 has 
been described in 'Ehr enreich et al.l (2011). In order to 
maximize the combination of wavelength coverage and 
sensitivity, near-UV spectra were observed with the STIS 
G230L/cenwave 2376 mode (At; « 600 km s~^) through 
the 52" X 0.1" slit. In order to constrain the effects of 
interstellar absorption on the stellar Mg II emission pro- 
file, we also acquired Mg II spectra with the STIS E230H 
mode {Av « 2.6 km s^^) through the 0.2" x 0.2" slit for 
GJ 832 and GJ 667C. 

Exposure times were generally short due to the pilot- 
study nature of this project; considerably longer inte- 
grations will be essential for future observations to mea- 
sure UV flare frequencies in exoplanet host stars and to 
produce high quality line profiles of important chromo- 
spheric and transition region tracers such as C II A1335, 
N V A1240, and Si IV A1400. A complete list of the M 
dwarf observations used in this work is presented in Ta- 
ble 1, and complete spectra for the four targets in which 
Lya, FUV, and NUV emission is observed are displayed 
in Figure 1. 

In order to place the M dwarf exoplanet host stars in 
context with other well-studied cool stars, we assembled 
archival spectra of AD Leo (M3.5Ve), HD 189733 (KIV), 
and the Sun (G2V). AD Leo is one of the best-studied 
fiare stars, and despite its extreme activity levels is has 
been the object of choice for models of extrasolar planets 
orbiting M dwarfs, mainly because of the lack of reason- 
able alternatives. We created a complete UV spectrum 
of AD Leo by combining STI S E140M monitoring ob- 
servations (jHawlev et al.ll2003l ) with an "average" NUV 



spectrum from lUE. The lUE spectrum is the average 
of flare and quiescent states, scaled to the FUV flux from 
STIS. HD189733 is an active K dwarf hosting the proto- 
typical transiting hot Jupiter. While high-quality COS 
spectra of HD 189733 exist in the archive (e.g., Linsky et 
al. 2012; Haswell et al. 2012), we elected to use the com- 
bined NUV-I-FUV (STIS E140M + E230M) spect ra of a 
surrogate star taken from the STARGAT library (jAvresI 
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Fig. 4. — Intrinsic, reconstructed Lya emission fluxes compared 
with the observed Mg II {h + k) flux, both scaled to the flux at 
0.16 AU from the star. The Lya and Mg II fluxes are observed to 
be well correlated. The slope of the relation is 14.8 ± 1.2. After 
correcting for the effects of interstellar absorption of the Mg II lines, 
the i<'(Lya)/F(Mg II) ratio for weakly active M dwarf exoplanet 
host stars is 10 ± 3. The expected Lyo emission line strength of 
GJ 1214 is shown as the open red circle. The observed GJ 1214 
Lya emission line flux upper limit is shown as the flUed red star. 



120101) to eliminate uncertainties about temporal variabil- 
ity between different observations. We chose HD37394 
(KIV variable star; V = 6.23) because it is a good match 
in stellar mass and activity level to HD 189733, and be- 
cause the FUV spectra of the two stars are qualitatively 
similar. The UV spectrum of HD37394 was then scaled 
to the V-magnitude of HD 189733 {V = 7.77). The 



spectrum of the quiet Sun was taken from iWoods et al.l 
([2009), and serves as the prototype main-sequence G- 
type star. 

4. INTRINSIC LYa PROFILE RECONSTRUCTION AND 
TIME VARIABILITY ANALYSIS 

4.1. Lya Reconstruction 

The stellar Lya emission line dominates the UV out- 
put of M dwarfs, containing approximately as much 
energy as the rest of the FUV-I-NUV spectrum com- 
bined ([France et al.l [2012af) . Therefore, measurements 
of Lya emission provide an important constraint on the 



6 



France et al. 



GJ832-temp 




-2x10"" [ , , I , , , \ iL^ , , , I , , ; 

2796.0 2796.5 2797.0 

Wavelength (A) 

Fig. 5.— STIS E230H profile of the Mg II k line 
(Ao = 2795.528 A) of GJ 832. The dotted magenta lines are the 
stellar emission and interstellar absorption components of the fit, 
respectively. Fitting the high-resolution echelle spectra of the Mg II 
k and h lines of GJ 832 and GJ 667 C, we estimate that the inter- 
stellar correction for the Mg II flux is ~ 30 — 35%. 

source term for calculations of the heati ng and pho- 
tochem istry of exoplanetary atmospheres (jLinskv et alJ 
l2012aD . Due to resonant scattering of neutral hydrogen 
and deuterium in the interstellar medium (ISM), the in- 
trinsic Lya radiation field cannot be directly measured, 
even in the nearest stars (logio A^(HI) < 18.5). In order 
to produce an accurate estimate of the local Lya flux 
incident on the planets around our target stars, the Lya 
profiles observed by STIS must be reconstituted. We 
do this using a technique that simultaneously fits the 
observed Lya emission line profile and the interstellar 
medium component, using the M PFIT routine to m ini- 
mize between the fit and data ()Markwardtll2009n . 

The intrinsic Lya emission line is approximated as a 
two-component Gaussian, as suggested for tr ansition re- 
gion emission lines from late-type stars by Wo od et alJ 
p997.) . This approximation is further justified by 
the reconstructed Lya profiles for M dwarfs, which 
show approximately Gaussian line shapes and little 
evidenc e for self-reversal compared to F, G, and K 
dwarfs (|Wood et al.|[2005l ). The Lya emission compo- 
nents are characterized by an amplitude, FWHM, and 
velocity centroid. Interstellar absorption is parameter- 
ized by the column density of neutral hydrogen (TV (HI)), 
Doppler b parameter, D/H ratio, and the velocity of the 
interstellar absorbers. We assumed a fixed D/H ratio 
(D/H = 1.5 X 10"^ Linsky et al. 1995; 2006). Degener- 
acy between the two emission components makes errors 
difficult to define for the individual fit parameters. We 
estimate the uncertainty on the integrated Lya fluxes 
from our technique by comparing the reconstructed Lya 
flux for a range of initial guesses on the parameters. The 
uncertainty on the integrated intrinsic flux is « 10 - 
20 % for the stars observed with E140M (GJ 832 and 
GJ 667C) and « 15 - 30 % for the stars observed with 
G140M (GJ 581, GJ 876, GJ 436). In order to test our 
methodology, we fitted the STIS E140M Lya profile of 
AU Mic. Our total integrated flux agrees with that pre- 
sented in Wood et al., (,2005i) to « 5%, with an identical 



derivation of the interstellar atomic hydrogen column on 
the line of sight (logio N(m) = 18.36 ± 0.01). 

Employing our iterative least-squares reconstruction 
technique, we also derived the intrinsic Lya spectrum 
of HD 189733. Using the 2010 Aprfl 06 STIS G140M 
spectra (L ecavelier d es Etangs et al, ,2012), we find an 
intrinsic fiux of F(Lya) = 7.5 x 10"^'^ erg cm~^ 
s^^ and an interstellar hydrogen column density logio 
iV(HI) 18.45. We also refitted the intrinsic Lya emis- 
sion from AD Leo, finding F(Lya) = 7.5 x 10^^^ erg 
cm^^ s^^ and an interstellar hydrogen column density 
logio N{m) = 18.47 ± 02. This value agr ees with the 
Lya reconstruction from iWood et all ()2005f) to w 30%, 
which gives us confidence in the iterative technique, even 
for sightlines with multiple interstellar velocity compo- 
nents. 

We present the individual Lya spectra in Fig- 
ure 2, and narrow/broad emission component ratios 
{Fnarrow / Fbroad) , spectral liuc widths, H I column 
densities, and Doppler 5- values in Table 2. We find 
Fnarrow/ Fbroad to be in the range 5-12 for the MUS- 
CLES stars, again independent of the resolution of the 
observations on which the reconstructions are based. For 
comparison, we find Fnarrow / Fbroad ^ 3 for the more 
active M dwarf (AD Leo) and more massive star (HD 
189733). Reconstructed Lya fluxes are given in Tables 
3, 4, and 5. GJ 1214 has no Lya emission detected, and 
we discuss this in greater detail in §5.1.1. 

4.2. Timing Analysis 

In our initial study of GJ 876, we described the 
large chromospheric/transition region flare observed in 
severa l UV emission line s during our COS observa- 
tions (jFrance et al.ll2012a[ ). We have carried out a simi- 
lar analysis on the entire MUSCLES sample to constrain 
the heretofore unknown level of UV variability in opti- 
cally quiet M dwarf exoplanet host stars. Light curves 
in several chromospheric and transition region lines were 
extracted from the calibrated three-dimensional data by 
exploiting the time-tag capability of the COS microchan- 
nel plate detector (France et al. 2010a). We extract a 
[Xi,yi,ti] photon list (where A is the wavelength of the 
photon, y is the cross-dispersion location, and t is the 
photon arrival time) from each exposure i and combine 
these to create a master [A,j/,t] photon list. The total 
number of counts in a [AA,Ay] box is integrated over a 
timestep At. We use a flux-dependent timestep of At 
= 40 ~ 120s for the MUSCLES targets. The instrument 
background level is computed in a similar manner, with 
the background integrated over the same wavelength in- 
terval as the emission lines, but physically offset below 
the science region in the cross-dispersion direction. The 
S/N was not high enough to temporally resolve the GJ 
1214 data. 

5. STELLAR UV EMISSION LINES AND THE FUV/NUV 
RATIO 

5.1. Chromospheric and Transition Region Emission 

Chromospheric and transition region emission lines are 
observed in all of the MUSCLES spectra, suggesting that 
many M dwarf exoplanet host stars have UV-active at- 
mospheres (Figure 1). The chromosphere is observed in 
H I Lya, C II AA1334,1335, Al II A1671, Fe II mufliplets 
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Fig. 6.— (top) HST-COS + STIS 1140 - 3140 A spectrum of 
GJ 1214. Emissions from C IV and Mg II are the only strongly 
detected features. At bottom, we show the ± 400 km s~^ spec- 
trum around H I Lya {upper panel, black histogram), C IV A1548 
{lower panel, red histogram), and N V A1239 {lower panel, blue his- 
togram) scaled up by a factor of 2 for display. The Mg II and C IV 
lines demonstrate the presence of an active atmosphere; several 
possibilities for the non-detection of Lya are presented in §5.1.1. 



near 2400 and 2600 A, and Mg II AA2796,2803; lines 
that are formed from ~ 1 - 3 x 10^ K. The transition re- 
gion is traced by lines with formation temperatures from 
~ 0.6 - 1.6 X 10^ K, including the C III 1175 multi- 
plet, Si III A1206, O V A1218, N V AA1239,1243, Si IV 
AA1394,1403, C IV AA1548,1550, and He II A1640. C II 
and C IV line fluxes for GJ 581, GJ 876, GJ 436 and GJ 
832 have been published in lLinskv et al.l (|2012aD . and a 
comprehensive set of stellar emission line measurements 
will be presented in a future work. In Table 5, we present 
the reconstructed Lya fluxes and the observed fluxes of 
C IV and Mg II. 



Figure 3 displays the C IV profiles for GJ 581, GJ 876, 
GJ 436, and GJ 832 scaled to the flux level at 0.16 AU 
from the parent star, a value near the center of the HZ for 
the MUSCLES targetfQ. While we conclude that most 
M dwarf exoplanet hosts have a non-zero local UV radi- 
ation field, the amplitude of this field is quite variable. 
The HZ C IV fiuxes from these four stars varies by a 
factor of 10 from the weakest (GJ 581) to the strongest 
(GJ 832) C IV emitter. This is particularly interest- 
ing in light of the established correlation between X-ray 
and C IV luminosity for M dwarf emission line (dMe) 
fiare stars (Bv rne fc Dovlelll98 9). With improved atmo- 
spheric models for M dwarfs or simultaneous UV and 
X-ray observations of a larger sample of M dwarf planet 
hosts, the C IV fiux in the HZ of low mass stars could be 
calibrated as a proxy for the EUV and sof t X-ray emis- 
sion in cident on the planetary atmospheres (jLinsky et al.l 
l2012a| ). C IV-based EUV irradiance estimates could pro- 
vide a complementary constr aint to those derived from 
coronal X-ray measurements (|Sanz- Forcada et aLll201lD . 

Lyo; is the brightest line in the UV spectrum of low- 
mass stars, although resonant scattering in the interstel- 
lar medium makes dir ect line profiles inaccess ible, even 
for the nearest stars (jWood et al.ll200(]l 120051 ). As de- 
scribed above, Lya reconstruction can be a powerful 
technique for determining the intrinsic stellar Lya line 
fiux, but this technique is not possible for more distant 
targets where the ISM completely removes stellar Lya 
photons from our line of sight (iV(Hl) > lO^'^ cm^^). In 
this case, one requires a proxy formed in a similar region 
of the chromosphere to estimate the local Lya emission 
flux. For our limited sample, we have found that Lya and 
Mg II are well correlated (see also Linsky et al. 2012a), 
and we present a scaling relation to infer the Lya fiux 
from a Mg II observation. 

Figure 4 shows the reconstructed Lya fiux versus the 
observed Mg II flux, scaled to 0.16 AU from the host 
star for GJ 581, GJ 876, GJ 436, GJ 832, and GJ 
667C (black filled circles). A linear fit to the correla- 
tion (i^(Lya) = a + bF{Mg II), where a = -2.3 ± 2.2 
and b = -1-14.8 ± 1.2 ) is also displayed in Figure 4. 
The observed F(Lya)/F(Mg II) ratio is the slope of this 
fit. Chromospheric Mg II emission lines of G and K 
dwarfs are known to show self-reversed profiles due to 
non-LTE effects in their upper atmospheres (e.g., Linsky 
& Wood 1996; Wood & Linsky 1998; Redfield & Linsky 
2002), howev er the situation is les s clear for M dwarfs. 
For example, I Redfield et al.l (|2002l ) present GHRS spec- 
tra of AU Mic, where no appreciable self-reversal of the 
reconstructed Mg II profile is observed. We therefore 
assume a simple Gaussian line shape for the M dwarf 
Mg II emission lines. Interstellar Mg II absorption will 
produce the appearance of a self-reversed emission line 
when the stellar velocity is close to the velocity of the 
Mg+-bearing interstellar cloud (e.g.. Wood et al. 2000), 
therefore a correction for interstellar Mg II absorption 
needs to be made in order to give an unbiased estimate 
of the F(Lya)/i^(Mg II) ratio for the MUSCLES stars. 

We cannot analyze the interstellar Mg II absorption 
in detail for most of our targets because the stellar and 



^ Ivon Braun et al.l S2m3 ] compute the HZ of GJ 436 to extend 
0.16 - 0.31 AU, so 0.16 AU may represent the inner edge of the HZ 
for this system. 
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interstellar features are unresolved in STIS G230L ob- 
servations. However, utilizing our STIS E230H high- 
resolution (Aw < 3 km s~^) spectra of GJ 832 and 
GJ 667C, the correction for interstellar Mg II absorp- 
tion can be estimated. The profile of the Mg II k line 
(Ao = 2795.528 A) of GJ 832 in shown in Figure 5. For 
GJ 832 and GJ 667C, we find that the interstellar ab- 
sorber removes 30 - 35% of the intrinsic Gaussian emis- 
sion line fiux. Therefore, we estimate that the scaling of 
the ISM-corrected Mg II emission line flux to the intrinsic 
Lya emission line fiux is « 10. Factoring in uncertainties 
on the interstellar correction and the Lya reconstruction, 
we estimate that this scaling relation is good to ~ 30%, 
or the ISM-corrected intrinsic J^(Lya)/J^(Mg II) ratio 
for weakly active M dwarf exoplanet host stars is 10 ± 3. 
Linsky et al. (2012a) find that the F(Lya)/F(Mg II) ra- 
tios for the more active M dwarfs AU Mic and AD Leo 
are 2.5 and 4.0, respectively. This suggests that the ratio 
may decrease with activity. 

5.1.1. Non-detection of Lya m GJ 1214 

GJ 1214 is the coolest (M6V, Te// 2949 ± 30 K; Kun- 
durthy et al. 2011) and most distant source in the MUS- 
CLES pilot study, and the only star where H I Lya is 
not detected. In the top panel of Figure 6, we display 
the full UV spectrum of GJ 1214. The only clearly 
detected stellar emission lines are those of C IV and 
Mg II. In the lower panel, we focus on the immediate 
spectral region around the lines of interest. The up- 
per limit from the coadded STIS G140M spectra of GJ 
1214 is F(Lya) <2.4 x lO^^^ erg cm^^ s^^. The short- 
wavelength component of the C IV doublet is shown 
in red, and we show a possible detection of N V in 
blue. Taking the observed Mg II flux from GJ 1214 
(2.21 (± 0.15) X 10~^^ erg cm~^ s~^), we can use the 
i^(Lya)/i^(Mg II) relation described above to compare 
the expected flux from Lya with the upper limit (red 
circle and red star in Figure 4). The upper limit on the 
GJ 1214 Lya flux is an order of magnitude below the 
predicted point based on the observed Mg II flux. Note 
however that this relation is for reconstructed Lya emis- 
sion lines and we have no information on which to base 
a reconstruction in GJ 1214. 

The presence of the UV metal line emission demon- 
strates that GJ 1214 maintains some basal UV fiux level. 
However, the non-detection of Lya is surprising in light 
of the results presented above. We suggest three possible 
explanations that might account for the non-detection. 
First, the 2012 COS and STIS G230L observations may 
have occurred during a period of stellar activity, whereas 
the star may have been UV-quict during the 2011 STIS 
G140M observations. This type of "binary" UV activity 
behavior ha s been observed pr eviously in the M8 fiare 
star VB 10 (iLinskv et al.llT995l) . Second, if the GJ 1214 
Lya emission line was significantly narrower than the 
Lya lines in the other M dwarfs, then it might be pos- 
sible that the somewhat longer interstellar sight line re- 
moved the signal. However, as shown by Wood et al.l 
(|2005D . neutral interstellar hydrogen columns rarely ex- 
ceed logio iV(HI) = 18.5 cm~^ for stars nearer than 50 
pc. Therefore, this second option is really that the stel- 
lar Lya emission line from GJ 1214 is unusually narrow. 
There is not enough information from C IV or Mg II to 
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Fig. 7. — (top) The total FUV (including the reconstructed Lya 
emission line) to NUV flux ratio is ~ 0.5 — 3 for M dwarfs and 
< 10^^ for K dwarfs and earlier. This ratio has important implica- 
tions for the oxygen photochemistry on Earth-like planets orbiting 
low-mass stars, {bottom) The fraction of the total UV energy in 
the Lyo line is > 35% for most of the M dwarfs (GJ 581, GJ 876, 
GJ 436, GJ 832, and AD Leo), < 8 % for GJ 1214, and < 1% for 
earlier- type stars. 



make a meaningful comparison of GJ 1214 line widths to 
the other stars in our sample. The third, more specula- 
tive, possibility is that for the low effective temperature 
of GJ 1214, a significant fraction of the available neu- 
tral hydrogen in the stellar atmosphere is locked up in 
molecules (mainly H2). As we will discuss in §5.3, we 
detect H2 emission in all of the other M dwarfs observed 
with COS. We do not detect H2 fluorescence in GJ 1214, 
but this could be the consequence of a feedback process 
driven by increasing molecular fraction. Lya is the en- 
ergy source for the observed fluorescence, so as fewer 
atoms are available to produce Lya, the H2 fiuorescence 
process will become unobservable in FUV observations 
even if more molecules are present in the stellar atmo- 
sphere. 

5.2. FUV/NUV Flux Ratios and the Lya Fraction 

In earlier work describing the UV spectrum of GJ 876, 
it was showi i that the FUV/NUV fiux ratio is > 10^ x the 
solar value (jFrance et al.ll2012al ). This is a consequence 
of an enhanced Lya line strength and a much lower NUV 
fiux owing to the cooler effective temperature of GJ 876. 
In Figure 7, we show that this is a generic property of M 
dwarfs. Here, we define the FUV fiux, F(FUV), as the 
total stellar flux integrated over the 1160- 1690 A band- 
pass, excluding Lya. Therefore the total FUV emission 
is SFUV = (F(Lya) + F(FUV)), the total FUV/NUV 
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ratio is 

T.FUV/NUV = {F{Lya) + F{FUV)) / F{NUV), (1) 

and the Lya fraction is defined as 

f{Lya) = F{Lya)/{F{Lya) + F{FUV) + F{NUV)), 

(2) 

where -F(Lya) is integrated over 1210 - 1222 A and 
F(NUV) is integrated over 2300 - 3050 A. 

The observed Lya, FUV, and NUV fluxes are pre- 
sented in Table 3. We put these values in context in 
Table 4 by presenting the Lya, FUV, and NUV lumi- 
nosities as the fraction of the bolometric luminosities. 
We calculate the bolometric luminosity for each star by 
creating a continuous spectrum from 0.115 - 2.5 Atm, us- 
ing t he appropriate M dwarf model atmosphere (jPicklesI 
|1998[ ) to cover the optical and IR portions of the SED. 
The model atmosphere flux is scaled to the STIS G230L 
data from 3050 - 3080 A. The use of the Pickles models 
will underestimate the bolometric luminosity somewhat 
(by not including the entire Rayleigh- Jeans tail into the 
mid-IR), however we expect that this effect is similar to 
or smaller than the uncertainty introduced by scaling to 
the NUV data, where the relative contribution of the 
photosphere is not clear in most cases. 

SFUV/NUV is in the range « 1 - 3 for aU M dwarfs, 
except GJ 1214. Mostly due to the lack of Lya emis- 
sion (§5.1.1), EFUV/NUV is in the range 0.50 - 0.62 
for GJ 1214. While the EFUV/NUV ratio is lower for 
GJ 1214 than for the rest of the M dwarfs in our sam- 
ple, it is still significantly higher than for earlier type 
stars. For HD 189733, this ratio has fallen to ^ 10~^, 
and it is < 10"'^ for the Sun. Similarly, /(Lya) spans 
37 - 75 % for most M dwarfs, is < 7.4% for GJ 1214, is 
« 0.7 % for HD 189733, and is w 0.04 % for the Sun. 
AD Leo has an absolute NUV fiux that is 4 - 40 times 
larger than for the less active M dwarfs in this survey, and 
the EFUV/NUV ratio could decrease dramatically dur- 
ing flares when stro ng Balmer continuum emiss ion domi- 
nates the spectrum (|Hawlev fc Pettersenlll991[ ). The im- 
plications of the large FUV/NUV flux ratio in M dwarfs 
are discussed in more detail below (§7.1). 

5.3. The Ubiquity of H2 Fluorescence: an Exoplanetary 

Origin ? 

We observe fluorescent H2 emission lines in GJ 581, 
GJ 876, GJ 436, and GJ 832. These features are pho- 
toexcited ("pumped") by Lya photo ns out of the v — 2 
level of the ground electronic state (jShullI IT978I ). Un- 
der the assumption of thermalized rovibrational popula- 
tions, this process is a signpost for 2000 - 4000 K molec- 
ular gas. We detect between four and eight fluorescent 
emission lines in each of these targets, pumped by Lya 
through the (1 - 2) R(6) A1215.73 A and (1 - 2) P(5) 
A1216.07 A transitions. The S/N of any of the indi- 
vidual lines is not sufficiently high for spectral fitting, 
so we coadded four to six of the brightest linetO- The 

The coadded H2 line profiles include at least four of the 
following emission lines: (1 - 7)R(6) A1500.24 A, (1 - 7)P(8) 
A1524.65 A, (1 - 6)R(3) A1431.01 A, (1 - 6)P(5) A1446.12 A,(l - 
7)R(3) A1489.57 A, and (1 - 7)P(5) A1504.76 A. 



coadded line profiles were normalized to unity and fit- 
ted by assuming a Gaussian emission line shape that 
has been convolved wi th th e 1500 A HST+COS line 
spread function (Kriss 20111) . The coadded fiuorescence 
lines and fits are shown in Figure 8. The H2 emission 
lines are consistent with the stellar radial velocities to 
within half of the 15 km s~^ wavelength solution un- 
certainty of the COS G160M mode, [vh2{581) = -5.5 
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cept in the case of GJ 581 (the lowest S/N H2 spectrum) 
the H2 emission features are spectrally unresolved. In 
this subsection, we describe three possible physical ori- 
gins for the observed fluorescence: 1) the stellar photo- 
spheres, 2) the atmospheres of the planets orbiting these 
stars, or 3) circumstellar gas disks populated by atmo- 
spheric mass-loss from the planets in these systems. 

Lya-pumped H2 emission lines w ere first detected i n 
sounding rocket spectra of sunspots (jJordan et al.lll977|) , 
and have been modeled assuming a thermalize d H2 pop- 
ulation at a temperature of 3200 K ()Shulllll978[ ). H2 fluo- 
rescence lines are also powerful diagnostics of the molec- 
ular surface layers of protoplanetary disks over a similar 
range of temperatures (Herczeg et al. 200j; [France et al.l 
2012b). These lines are not clearly detected in previ- 
ous STIS observations of M dwarfs without planets (AD 
Leo, Proxima Gen, EV Lac, and AU MicEl), although 
this could be related to much larger instrumental back- 
grounds associated with the STIS echelle modes. In any 
event, the 3200 K sunspot temperature is characteristic 
of the photospheric effective temperature of the MUS- 
GLES stars, and simple energy budget calculations (e.g., 
France et al. 2012a) suggest that the fluorescence may 
be consistent with an origin in a hot molecular stellar 
surface layer. 

However, it is interesting that these lines are ei- 
ther not detected or are marginally detected with much 
lower H2/atomic line ratios in the non-planet hosting 
stars. H2 is the primary constituent of gas giant plan- 
ets. After H I Lya, the electronically excited emission 
spectrum of H2 do minates the FUV spectrum of giant 
planets jBr oadfoot et al.lll979l:rWolven fc Feldmanlll998l: 
iGustin e~al..,2010,) . Fluorescence driven by solar Ly/3 is 
the main constituent of the daygl ow (non-auroral illu- 
minated disk) spectrum of Jupiter (jFeldman et al.|[l993f) 
and in areas where the Jovian surface temperature is el- 
evated (e.g., the Shocmakcr-Levy 9 impact site; Wolven 
et al. 1997), Lya-pumped fluorescence becomes strong. 

All four of the stars showing strong H2 fluorescence 
host planets of Neptune to super- Jupiter mass. Our anal- 
ysis cannot rule out the possibility that the ubiquitous 
Lya-pumped H2 fluorescence observed towards M dwarf 
planet hosts is the result of photoexcited dayglow emis- 
sion in planetary atmospheres. In light of the recent ten- 
tative detections of FUV H2 emissio n from the plan etary 
systems around HD 209458 ( Fran ce et all l2010bD and 



brown dwarf 2MASS J 12073346-3932539 ( [France et all 



AU Mic displays H2 fluorescence from lower temperature 
molecular gas, but this g as is consistent with an origin in a remnant 
protoplanetary gas disk ll'Vance et al.|[2007 V 
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TABLE 3 

M DWARF Broadband UV and Reconstructed Lyq Fluxes (ergs cm~^ s"-')''. 



Target 


d (pc) 


F(Lya)b 


F(FUV)'= 


F(NUV)'^ 






(AA = 1210 — 1222 A) 


(AA = 1160 — 1690 A) 


(AA = 2300 - 3050 A) 


GJ 581 


6.3 


3.0 X 10"" 


3.6 X 10-'^^ 


2.3 X 10-13 


GJ 876 


4.7 


4.4 X 10"" 


1.7 X 10-13 


3.9 X 10-13 


GJ 436 


10.2 


3.5 X 10"" 


3.3 X 10-14 


2.8 X 10-13 


GJ 832 


4.9 


5.0 X 10~^2 


9.8 X 10-1" 


2.2 X 10-12 


GJ 667G 


6.9 


7.6 X 10"" 






GJ 1214 


13.0 


< 2.4 X 10-15 


1.0 X 10-14 


2.0 X 10-1" 


AD Leo 


4.7 


7.5 X 10"^^ 


3.4 X 10-^^ 


9.3 X 10"^^ 


HD 189733° 


19.5 


7.5 X 10"" 


4.1 X 10-13 


1.0 X 10^1° 


Sun 


1 AU 


5.9 X 10° 


5.9 X 10° 


1.6 X 10"' 



^ Flux measurements are averaged over all exposure times. 

^ Uncertainty on the reconstructed Lya flux is estimated to be between 10 - 20% for the stars observed with E140M 
(GJ 832, GJ 667C, and AD Leo) and 15 - 30% for stars observed with G14GM (GJ 581, GJ 876, GJ 436, and HD 
189733). 

" F{F\JV) refers to the integrated 1160 - 1690 A flux excluding Lya. Flux uncertainties for the broad band measure- 
ments are dominated by instrumental calibrations and target acquisition errors for faint sources. The flux errors for 
i^(FUV) and F(NUV) are estimated to be < 10% and < 5%, respectively. 
'I AA is the bandpass over which the flux is integrated. 

HD 189733 Lya reconstruction made from direct STIS G140M observations. F(FUV) and i^(NUV) were measured 
using a proxy star (HD 37394) scaled to the V-magnitude of HD 189733. 



TABLE 4 

M DWARF Broadband UV Luminosity (ergs s-i)^ 



Target 


d (pc) 


logio L(Lya) 


logio L(FUV)*= 


logio L(NUV) 






(AA'^ = 1210 - 1222 A) 


(AA = 1160 - 1690 A) 


(AA = 2300 - 3050 A) 


GJ 581 


6.3 


27.16 


26.23 


27.03 


GJ 876 


4.7 


27.07 


26.65 


27.01 


GJ 436 


10.2 


27.65 


26.62 


27.55 


GJ 832 


4.9 


28.15 


26.44 


27.80 


GJ 1214 


13.0 


<25.69 


26.31 


26.61 




logio L{Lya) / Ltoi'^ 


logio L(FUV)/L6„i 


logio L{^Vy)/Ltoi 


GJ 581 


6.3 


-4.11 


-5.04 


-4.24 


GJ 876 


4.7 


-4.26 


-4.68 


-4.32 


GJ 436 


10.2 


-4.07 


-5.10 


-4.17 


GJ 832 


4.9 


-3.50 


-5.21 


-3.85 


GJ 1214 


13.0 


<-6.24 


-5.62 


-5.31 



Flux measurements are averaged over all exposure times. 
^ X(FUV) refers to the integrated 1160 - 1690 A flux excluding Lya. 
" AA is the bandpass over which the flux is integrated. 
=1 Lhoi is the integral of the X-ray through IR flux of each star (§5.2). 



I2010al) . it is tempting to attribute an exoplanetary origin 
to these lines. We caution that further study is required 
to differentiate between these two scenarios. If the H2 
in the GJ 436 system originates in the atmosphere of 
GJ 436b, we would have expected to see a significant ve- 
locity shift of the fluorescent lines relative to the stellar 
radial velocity (■i;436b = 110 - 120 km s~^ at orbital phase 
« 0.8). However, the H2 lines are consistent with the ra- 
dial velocity of the star, indicating little direct dayglow 
and/or auroral contribution from GJ 436b. The upper 
limit on the H2 (1 - 7)R(3) emission line flux at -f 110 km 
(assuming unresolved emission lines) is « 6 x 10^^*" 
erg cm~^ s~^ in the GJ 436 spectrum. Modeling of both 
the H2 fluorescence spectrum from M dwarf photospheres 
and the expected UV emission spectrum from nearby M 
dwarf planetary systems would be extremely valuable. 
A final possibility is that the observed H2 fluorescence 



originates in a circumstellar gas envelope that is being re- 
plenished by atmospheric mass-loss from escaping plan- 
etary atmospheres. The narrow H2 line-width in GJ 
436 rules out the possibility that the fluorescing gas is 
confined to the semi-major axis of the planetary orbit 
(a = 0.03 AU for GJ 436b), however over the - 5 Gyr 
lifetimes of these systems, this material could disperse 
into a low density "second generation" circumstellar disk. 
While detailed modeling of such a disk is beyond the 
scope of this work, we can make a rough estimate of the 
planetary mass-loss rate required to sustain the observed 
H2 emission. 

The balance equation of H2 photodissociation and re- 
plenishment is (3dissn{il2) = ^loss- Due to the negli- 
gible 912 - 1120 A continuum flux from inactive M- 
type stars, we assume that the interstellar radiation field 
(ISRF) dominates H2 dissociation. The photodissoci- 
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Fig. 8. — Coadded H2 emission spectra detected in the MUSCLES pilot program. Each velocity profile is the coaddition of four to 
six fluorescent emission lines. The dotted line is the intrinsic Gaussian line profile which has been convolved with the COS line spread 
function (solid blue line) for comparison with the data (black histogram). The red error bars are representative of the uncertainty on the 
normalized line profiles. The origin of these lines is not yet clear; more detailed modeling is required to differentiate between a photospheric 
and exoplanetary origin for this hot, photoexcited molecular gas. 



ation rate is (3diss = BdissPo- Bdiss is the dissocia- 
tion fraction of molecules excited into the Lyman and 
Werner bands following absorption of a 912 - 1300 A 
photon, K, 0.2 for a hot H2 gas. /3o is the standard H2 
photoabsor ption rate from the ISRF, /3o = 5 x 10"^" 
s-i (iJm^ Il975h . We assume an H2 column den- 
sity of Af(H2) « 10^^ cm-2 (jFrance et all l2012ar) and 
a fiducial cloud diameter of dH2 = 1 AU such that 
n(H2) = N{Yi2)ldH2 « 7 X 10^ cm'^. The replenish- 
ment term from the escaping atmospheres of our short- 
period planets is $;oss = Mioss iri~^l F(H2/HI) V^^f, 
where Mioss is the planetary atmosphere mass-loss rate, 
i^(H2/HI) is the fraction of the total cloud in molecu- 
lar form, and Vemit is the volume of the emitting area 
iVemit — '^rjj2ZH2)- Modcls of cscapiug hot Jupiter 
atmospheres find that most of the escaping material is 
atomic (e.g., Tremblin & Chiang 2012), so we assume a 
low molecular fraction F(H2/HI) ~ 10^^. Any disk pop- 
ulated by an orbiting planet will most likely be geomet- 
rically thin, and we take zh2 = 0.01 dH2 as the vertical 
height. Keeping the crudeness of this approach in mind, 
we find that the H2 disk can be sustained for mass-loss 
rates Mioss ^ 6 X 10^" g s~^, well within the range of 
mass- loss estimat es of short-period planets around G, K, 
and M-tvpe stars (iVidal-Madiar et al.ll2008l:lLinskv et ah 
l2010l : lLecavelier des Etangs et al.ll2012l : lEhrenreich et al. 



l2011l:lEhrenreich fc Desertll2011[ ). Therefore, the H2 fluo- 
rescence observed in our MUSCLES sample appears com- 
patible with either a photospheric or circumstellar origin. 

6. TIME VARIABILITY OF M DWARF UV RADIATION 
FIELDS 

We show examples of the UV light curves of GJ 581, 
GJ 876, GJ 436, and GJ 832 in Figure 9. The gaps 
in the light curves shown in Figure 9 are due to Earth 
occultations during the orbit of HST. We observe time- 
variability in the UV spectra of all of our M dwarf ex- 
oplanet host stars observed with COS (GJ 1214 has in- 
sufficient S/N for timing analysis), ev en on weakly ac- 
tive stars (IW alkowicz fc Hawlevl[2009l ). In Fi gure 9, we 
show example light curves from all four targets in chro- 
mospheric and transition region tracers Si III A1206, N V 
AA1239,1243, C II AA1334,1335, Si IV AA1394,1403, and 
C IV AA1548,1550, and the detector background, evalu- 
ated in At = 40s or 120s timesteps. Our observations 
were designed to provide a first estimate of the flux lev- 
els and temporal variability in low-mass exoplanet host 
stars, so the observations do not span a sufficient baseline 
to determine the UV flare frequency. However, we can 
qualitatively identify three fiare behaviors in our spectra: 

1 ) The Impulsive Flare - The GJ 876 and GJ 436 UV 
flares exhibit a qualitatively similar behavior to optical 
light flares on M dwarfs (e.g., Figure 5 of Kowalski et 
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Fig. 9. — We have detected a range of temporal variability scales in the chromospheric and transition region emission lines from weakly 
active M dwarf exoplanet host stars. We separate the time variability into three qualitative groups: (i) The Impulsive Flare, shown in 
several chromospheric and transition region tracers (C II, Si III, C IV, and N V) in GJ 876 and GJ 436 (top row), (ii) The Slow Rise 
Flare, shown in Si IV in GJ 832 (lower left), {iii) The Short Flare, shown in C IV in GJ 581. This demonstrates that time- variable, 
active atmospheres are a generic property of M dwarf exoplanet host stars. The gray vertical dotted lines denote the end of an individual 
H ST-COS exposure. Exposure IDs are included. 



al. 2009). We observed a flare peak/quiescent flux ratio 
of > 10 on GJ 876, with an approximately exponential 
decay (rgre ^ 400s). The UV flare behavior on GJ 876 is 
similar t o the factor of ^ 5 - 20 flux increase observed on 
AD Leo (jHawle v ct al."2003'). GJ 876 also displays X-ray 
flare activity (Poppcnhacgcr et al. 2010). A similar but 
lower amplitude C IV flare was observed on GJ 436. The 
C IV flux approximately doubled in 240s, then decayed 
to pre-flare levels on a timescale r436 ~ 1200s. 

2) The Slow Rise Flare - The Si IV light curve of GJ 
832, binned to At = 120s timesteps, is shown in the lower 
left hand panel of Figure 9. Compared with the abrupt 
flux increase seen in GJ 876 and GJ 436, the GJ 832 Si IV 
flux increases by a factor of 1.5 - 2.0 over the course of 
the first 1200s exposure. The peak and turnover are lost 
to Earth occultation between COS exposures for GJ 832, 
but the flare is observed to slowly decay during the last 
960s of G130M observation. 

3) The Short Flare - GJ 581 displays a shorter duration 
flare event, shown in the lower right panel in Figure 9. 
The C IV light curve shows a flux increase of a factor 
of ~ 2 - 3 at Texp ~ 750s. The signal-to-noise in the 
GJ 581 UV emission lines is low, but one can resolve 
the flare using At = 40s timesteps. This 200s long C IV 
emission enhancement comes during the first HST orbit, 
when the basal transition region flux was approximately 



a factor of two larger than during the second HST orbit 
following Earth occultation. 

The conclusion from this section is that even with a 
small sample of relatively low S /N spectra obtained over 
limited temporal baselines, we have identified variability 
with a range of amplitudes and durations. Continuous 
HST observing campaigns, spanning 6-10 orbits each, 
of a larger sample of exoplanet host stars will be essen- 
tial to better constraining the frequency, duration, and 
amplitude of the variable UV radiation environment of 
planets around M dwarfs. Alternatively, long-duration 
NUV observations could be made with balloon-borne ex- 
periments. Time-domain studies at optical wavelengths 
are also of interest for characterizing the flare-frequency 
of the largest flares that are panchromatic events (e.g., 
Hawley et al. 2003). This is particularly relevant in 
light of the projected paucity of UV-capable observato- 
ries. A more complete understanding of UV variability 
on weakly active and inactive M dwarfs, in combination 
with better constrained flare frequencies for stars of a 
variety of activity strengths, would allow one to use the 
wealth of data from current and upcoming large optical 
surveys to inform planetary atmosphere models. 



7. DISCUSSION 
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7.1. UV Radiation Fields for Atmospheric 
Photochemistry Models 

Owing to the scarcity of M dwarf exoplanet host spec- 
tra previously available in the literature, authors mod- 
eling the atmospheres of Earth-like planets around M 
dwarfs have typically been drive n to one of two extremes: 
assuming no UV flux ( Segura et al.l l2005t iRauer et al.l 
120111 : iKaltenegger et all i2011i) or assuming the spec- 
trum of one of the mos t active M dwarfs known , AD 
Leo (jSegura et al.l I2Q05I : iWordsworth et al.l |2010[) . In 
the case of no UV flux, the short wavelength spectral 
cut off is determined by the photosphere of a theoreti- 
cal star with Teff < 3500 K, typically having negligible 
flux b elow 2500 A and no chr omospheric emission fea- 
tures. iWalkowicz et al.l ()2008() demonstrated that even 
M dwarfs with modest X-ray emission and no chromo- 
spheric Ha emission display NUV spectra that are qual- 
itatively similar to more optically active stars. Follow- 
ing this work, Kaltcnegger et al] (1201 ID noted that an 
active UV atmosphere is most likely present on almost 
all M dwarfs, which is consistent with the observations 
presented here. Additional observations are required to 
understand the UV flux from M dwarfs as a function of 
stellar age and metallicity, for spectral types later than 
M6, and for a larger sample of M dwarfs without X-ray 
emission or Ha activity. 

The MUSCLES observations show that the UV spec- 
tral behavior of four of our stars (GJ 581, GJ 876, GJ 436, 
and GJ 832) is similar to AD Leo in both the spectral 
and the temporal sense. The MUSCLES stars (except 
GJ 1214) showed similar a SFUV/NUV flux ratio and 
/(Lya) to that of AD Leo, and they also displayed time 
variability with UV emission line amplitudes changing by 
factors of > 2 during our brief observing snapshots. 

While AD Leo appears to be a more appropriate choice 
for exoplanetary modeling input, a primary difference is 
the absolute flux of the NUV radiation field from AD 
Leo. AD Leo has factors of ~ 4 - 40 stronger abso- 
lute flux levels at its effective HZ distance than do the 
rest of the MUSCLES targets. Furthermore, the strong 
Balmer continuum associated with its frequent flare ac- 
tivity will periodically drive the SFUV/NUV ratio sig- 
nificantly below the level (^1-3) found for the less ac- 
tive M dwarfs. This could lead to significant differences 
between the atmospheric chemistry of a planet orbiting 
an AD Leo-like star and one orbiting a more typical M 
dwarf. It has been recently shown that the strong FUV 
flux and high SFUV/NUV flux ratio incident on a ter- 
restrial planet orbiting in the HZ of an M dwarf can lead 
to the buildup of a large, abiotic atmospheri c O2 abun- 
danc e through the photodissociation of CO2 ()Tian et al.l 
I2012D . The creation of O3 via O2 photolysis depends sen- 
sitively on the strength of the 1700 - 2400 A radiation 
field and the balance of the ocean-atmosphere system on 
these planets. In some models, this process can lead to 
detectable levels of O3 without the presen ce of an active 
biosphere (jDomagal- Goldman et al.ll2012f) . 

Spectral observations of O2 and O3, along with CH4 
and CO2, are expected to be the most important signa- 
tures of biol ogical activity on planets with Earth-like at - 
mospheres (|Des Marais et al.l 120021: iSeager et al.ll2009D . 
which may be the case for such planets orbiting solar- 
type stars. However, the above work has shown that 
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Fig. 10. — Photo-absorption cross-sections of H2O (green dashed 
line; Mota et al. 2005) and CO2 (purple dashed line; Yoshino et al. 
1996) showing the spectral overlap with chromospheric and tran- 
sition region emission lines from an M dwarf (black histogram). 
Bright chromospheric and transition region emission lines are la- 
beled. 

when realistic inputs for the UV radiation field are in- 
cluded, it may be possible to produce O2 and O3 with 
observable abundances on terrestrial planets in the HZ of 
M dwarfs without the presence of biological life. There- 
fore, a detailed knowledge of the stellar spectrum will 
be critical to interpreting the observations of potential 
biomarkers on these worlds when they are detected in 
the coming decades. Until such time that stellar models 
of M-type stars can reliably predict the observed UV- 
thru-IR spectrum, a direct UV observation will be the 
most reliable means for determining the radiation field 
incident on planets orbiting M dwarfs. 

7.2. UV Transit Spectroscopy of M dwarf Planets 

In Section 6, we described the range of time- variability 
observed in the UV chromospheric and transition region 
emission lines from our M dwarf exoplanet host stars. 
We observe emission line flux variability of 50 - 500% 
on timescales shorter than the typical transit duration 
of short-period exoplanets (few hours). For comparison, 
the typical depth of a resolved UV transit observation is 
~ 5 - 15% (Vidal-Madiar et al. 2003: Bcn-Jaffcl 20071; 
iLinskv et al.l 120101 : iLecaveher des EFangs et al.l L2012). 
The importance of stellar variability for UV transit 
studies of G and K dwa rfs is just now being ex- 
plored ('Haswe ll et al.ll2"012D . While an older solar- type 
transiting planet host star such as HD 209458 has been 
shown to have very steady transition region emission 
(observed in Si IV; Linsky et al. 2010), UV variabil- 
ity on younger, less-massive transiting planet hosts like 
HD 189733 may significantly hinder the ability to detect 
transits and characterize their depth and ingress/egress 
profiles. GJ 436 and GJ 1214 both host transiting planets 
(Neptune and super-Earth mass planets, respectively), 
and our results suggest that transit studies using UV 
chromospheric emission lines as background sources to 
search for signs of atmospheric escape may be compli- 
cated by stellar activity. NUV searches for O3 in the 
atmospheres of Earth-like planets will also likely be im- 
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pacted by M dwarf stellar activity, even for weakly active 
stars. Long-baseline observations over multiple transits 
will be required to mitigate the effects of stellar activity 
in future M dwarf exoplanet transit observations. We 
will consider this issue in more detail in a future work. 

7.3. The MUSCLES Spectral Database 

As part of the MUSCLES pilot program, we have 
started building a library of medium/high resolution UV 
spectra of M dwarf exoplanet host stars. These data are 
available to download for use in photochemical models of 
Earth-like, ice giant, and gas giant planetj*^. 

The high-resolution of the MUSCLES sample is impor- 
tant for detailed molecular emission/absorption model- 
ing. As most of the UV flux from low-mass stars arises 
in discrete emission lines, molecular species in the atmo- 
spheres of HZ planets will be "selectively pumped" : only 
species that have spectral coincidences with stellar emis- 
sion lines will be subject to large energy input from the 
host star. H2O, CH4, and CO2 in particular are highly 
sensitive to far-UV radiation, and provide an illustration 
of the utility of spectrally resolved data. CH4 absorbs 
in a broad band at A < 1500 A, and the photoabsorp- 
tion rates are relatively insensitive to the spectral res- 
olution. CO2 displays a banded structure over a broad 
absorption region from 1200 - 1600 A, and high resolu- 
tion data show a spectral coincidence between the narrow 
CO2 bands and the Si III A1298 + O I A1304 muhiplet 
and the C II A1335 stellar emission lines. H2O dissoci- 
ation is influenced by far-UV photons as the molecules 
absorb in a broad continuum from 1400 - 1900 A, B - X 
absorption bands near 1280 A, and sharp Rydberg tran- 
sitions at A < 1240 A t hat overlap with the Lya emis- 
sion line (|Yi et al.l 120071 ). The spectral overlap between 
the photo-absorption cross-sections of CO2 and H2O and 
the chromospheric and transition region emission lines 
from an M dwarf are shown in Figure 10. Low-resolution 
spectra or molecular absorption cross-sections will fail to 
preserve the specific energy absorption in these narrow 
(AA < 10 A) wavelength coincidences and could lead to 
an underestimation of the dissociation rates of important 
atmospheric species. 

8. SUMMARY 

We have presented the first study of the UV radiation 
fields around nearby M dwarf planet hosts that covers 
both FUV and NUV wavelengths. Using a combination 
of data from the COS and STIS instruments on HST, 
we present 1150 - 3100 A spectra for GJ 581, GJ 876, GJ 
436, GJ 832, GJ 1214, and Lya and Mg II observations 
for GJ 667C. Using an iterative least-squares approach, 
we reconstruct the intrinsic Lya radiation field strengths 
assuming a two-component Gaussian emission line. Us- 
ing the time-tagged capability of the COS FUV detector, 
we calculated emission line light curves in several chro- 
mospheric and transition region diagnostics. The main 
results of this work are: 

1. All six of the stars in our sample display both FUV 
and NUV emission line flux, even though these 

12 The combined 1150 - 3100 A spectra (0.05 A pixel"! 
binning) and the full-resolution, coadded COS spectra 
of the MU SCLES targets are availab le for download at: 
|http : //cos . Colorado . 6du/~kevinf /muscles .html | 



stars do not display Ha emission (however, they 
do show some degree of Ha absorption and Ca II 
emission lines) . This suggests that an active stellar 
UV atmosphere should be included in calculations 
of atmospheric chemistry on planets orbiting most 
M dwarfs. 

2. The Lya emission lines are well-flt by a two- 
component Gaussian emission line model with nar- 
row (FWRMnarrow ^ 90- 130 km s-1) and broad 
(FWUMbroad ^ 160 - 300 km s'^) components. 
The Fnarrow I Fbroad Aux ratio raugcs from ~ 5 - 
12. Interstellar H I column densities are derived; 
18.06 < logioA^(HI) < 18.47 cm^^ Jqj. q^j^. targets 
{d = 4.7 - 10.3 pc). 

3. Lya emission was not detected towards the M6 
super-Earth host star GJ 1214, although emission 
from Alg II and C IV indicate that an active atmo- 
sphere is present, at least during certain times. 

4. The reconstructed Lya line flux correlates with the 
observed Mg II line flux. The ISM-corrected in- 
trinsic F(Lya)/F(Mg II) ratio for weakly active M 
dwarf exoplanet host stars is 10 ± 3. 

5. The total FUV/NUV flux ratio is in the 
range w 0.5 - 3 for all M dwarfs, including AD Leo. 
This ratio is < lO"'^ for KIV dwarfs and < 10"'^ for 
the Sun. The percentage of total UV energy from 
the star in the Lya emission line is: 1) 37 - 75 % for 
GJ 581, GJ 876, GJ 436, GJ 832, and AD Leo, 2) < 
7.4% for GJ 1214, 3) « 0.7 % for KIV dwarfs, and 
4) K. 0.04 % for the Sun. The FUV/NUV flux ra- 
tio can have a profound effect on the production of 
O2 and O3 in the atmospheres of Earth-like planets 
orbiting M dwarfs. 

6. Fluorescent emission from hot H2 (T(H2) ~ 2000 - 
4000 K) is detected towards all four stars with mod- 
erate S/N COS spectra (GJ 581, GJ 876, GJ 436, 
and GJ 832). A flrst order calculation suggests that 
this emission originates in the stellar photosphere, 
but this emission could also arise on the illuminated 
dayside of the planets in these systems or in a cir- 
cumstellar envelope replenished by mass-loss from 
the planetary atmosphere. 

7. Time-variability is detected in the chromospheric 
and transition region emission lines observed from 
all four stars with moderate S/N in the COS spec- 
tra (GJ 581, GJ 876, GJ 436, and GJ 832). Vari- 
ability with amplitudes 50 - 500% is observed on 
time scales of 10^ - 10'^ s. This will most likely 
make UV transit studies of these systems challeng- 
ing. 
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